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Abstract. Systems of globally coupled logistic maps (GCLM) can display complex collective behaviour 
characterized by the formation of synchronous clusters. In the dynamical clustering regime, such systems 
possess a large number of coexisting attractors and might be viewed as dynamical glasses. Glass properties 
of GCLM in the thermodynamical limit of large system sizes A*' are investigated. Replicas, representing 
orbits that start from various initial conditions, are introduced and distributions of their overlaps are 
numerically determined. We show that for fixed-field ensembles of initial conditions all attractors of the 
system become identical in the thermodynamical limit up to variations of order l/\/iV, and thus replica 
symmetry is recovered for TV ^ oo. In contrast to this, when random-field ensembles of initial conditions 
are chosen, replica symmetry remains broken in the thermodynamical limit. 

PACS. PACS-05.45.-aNonlinear dynamics and nonlinear dynamical systems - PACS-05.45.Xt Synchro- 
nization; coupled oscillators - PACS-75.10.Nr Spin-glass and other random models 

1 Introduction tiation |^. Understanding the properties of GCLM can 

be seen as a first step towards grasping the dynamics and 
The rich collective behaviour displayed by globally cou- emergent properties of real, high-dimensional systems. The 
pled logistic maps (GCLM) (l[| has made them to be- dynamical equations describing the system are 
come a paradigm of complex dynamical systems. Initially, 

GCLM were introduced as a mean field approach to cou- n 

a:.(t+l) = (l-e)/(x,(t)) + ^^/(x,(t)) (1) 
pled map lattices. Subsequently, they have been used as -"^ j^i 

metaphors of neural dynamics, ecology, and cell difFeren- where 
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systems, where the presence of quenched disorder causes 
f(x) = 1 — ax^ (2) frustration and a large number of macroscopic configura- 

tions are possible 0. For this reason, GCLM have been de- 
is the logistic map. The parameter e e [0, 1] gives the 

scribed as a dynamical counterpart to spin glasses Q , [|0[ . 
strength of coupling among elements. For e = the ele- 

In a previous publication [|l2| two of us have introduced 
ments evolve independently, and for e = 1 they are syn- 

a replica description for this system, defined overlaps and 
chronized already after the first iteration and follow iden- 

numerically tested replica-symmetry breaking and ultra- 
tical trajectories ever after. Between these two extreme be- 

metric properties of GCLM. Our analysis has revealed a 
haviours, a broad spectrum of collective dynamics emerges. 

strong size dependence of collective dynamics, indicating 
The dynamics is strongly sensitive to the control parame- 

that replica symmetry might be recovered in the thermo- 
ter a of the logistic map and depends on the size N of 

dynamic limit A^ — > oo. The aim of the present paper 
the system. Figure 1 shows a rough phase diagram of 

is to investigate systematically the asymptotic statistical 
GCLM, based on the collective behaviour of the system 

properties of GCLM in this limit, 
which is reached after (sometimes, very long) transients 



||l|,y|5|]. The diagram includes both the parameter interval Our main result is that the asymptotic behaviour ob- 
a < 1.4 where dynamics of an individual map is peri- served in the thermodynamic limit is strongly dependent 
odic and the interval 1.4 < a < 2 with chaotic individual on how the ensemble of initial conditions is prepared. In 
dynamics. It contains two large domains of synchronous previous studies |l^,|l^, the procedures used for random 
and turbulent phases. They are separated by a region with generation of initial conditions had a special property: in 
glass-like behaviour. In the synchronous domain the states the limit A^ ^ oo all generated initial conditions were 
of all elements are identical and the ensemble has the same effectively identical up to order 1/\/N. Therefore, all ex- 
dynamics as a single map. In the turbulent phase, the en- plored attractors in the glass phase became equivalent up 
semble of maps is essentially desynchronized, though non- to variations of order 1/Vn and the replica symmetry 
trivial correlations have been detected even there ||[. The was recovered in the thermodynamic limit. Now we show 
glass region is characterized by the formation of various that, if the initial conditions are prepared in such a way 
dynamical clusters. that the initial field always retains macroscopic fluctua- 
Numerical simulations of GCLM have shown that, in tions, the replica symmetry is clearly broken in the ther- 
the dynamical glass phase, the system displays sensitivity modynamic limit N ^ oo. Thus, GCLM represent the 
to initial conditions. For fixed parameters e and a (and first known example of a dynamical glass with replica- 
given N), a multiplicity of attractors can be reached [|l], symmetry breaking and, as we show, this important sta- 
|7|. This property is similar to what is observed in glassy tistical property does not represent a finite-size effect. Our 
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2 Characterization of attractors 
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Fig. 1. Rough phase space of GCLM showing the main three 
phases of the system. 

results also imply that the boundaries between different 
regions illustrated in Fig. 1 depend on the ensemble of 
initial conditions, and that, for broad ensembles and e not 
too large, fully synchronized attractors can coexist with 
multi-cluster attractors. 

In the next section we introduce dynamical and statis- 
tical measures needed to characterize clustered states of 
GCLM. They include the splitting exponent, earlier pro- 
posed by Kaneko |ll[ , and a new repulsion exponent which 
we suggest. Replicas and their overlaps are defined and 
attraction basin weights are considered here. In Section 3 
we perform a detailed analysis of the role of initial condi- 
tions. Replica symmetry breaking and ultrametric prop- 
erties of GCLM in the thermodynamic limit are investi- 
gated in Section 4 under a truly random choice of initial 
conditions. The paper ends with a discussion of the main 
results, which are compared to the properties of other dis- 
ordered systems. 



An attractor of the dynamical system ffl) is characterized 
by the formation of a number of clusters out of the ini- 
tially symmetrical ensemble of maps. This is one of the 
most intriguing properties of GCLM. Within each cluster 
all elements are completely synchronized. A partition of 
N maps into /C clusters is defined by indicating the num- 
bers Nk of elements in each cluster, k — 1, . . ./C. In the 
following, we assume that the clusters have been ordered 
such that iVi > 7V2 > . . . > Nk:-i > N^.. Even if only 
two clusters are present, this can correspond for N —f oo 
to a huge number of different partitions, since the relative 
sizes of clusters may vary. 

In the periodic region, and for e = 0, the elements can 
be trivially classified into N/P groups, where P is the pe- 
riod of the single map, and elements within each group 
follow the periodic orbit of the single map with different 
phases. At large enough e the number of simultaneously 
stable clusters decreases and their dynamical behavior dif- 
ferentiates. This happens approximately in the whole area 
labeled "dynamical glass" in Fig. 1. For e large enough, all 
of the maps are synchronized and the dynamics reduces to 
that of the single element (synchronous phase). Note that 
fully synchronized attractors and multi-cluster attractors 
can coexist in some region of parameter space. 

In the glass phase, attractors of GCLM correspond 
to different partitions in clusters. The global dynamics of 
each attractor can be periodic, quasiperiodic or chaotic. 
The periodic collective dynamics is by far the most com- 
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mon, and is typically found even in the parameter region "clusters" of a single element) in the chaotic domain of 

where the dynamics of a single map is chaotic. the logistic map. 

Due to the unavoidable finite precision of digital com- 

2.1 Splitting and repulsion exponents puters, the simulation of the deterministic system p) can 

lead to pseudo-attractors which are not stable against 
The route to synchronization can be easily found by study- 
transversal perturbations. In the results to be presented, 

ing how the distance between elements evolves in time. 

we have computed Af, for all orbits and discarded unstable 
This is ruled by the equation 

attractors. 

There is another condition which must be required for 

Xi{t + 1) — Xj{t + 1) = — (1 — e)a (xf{t) ~ x'^it)) ■ (3) having a stable partition and which, to our knowledge, has 

not been made explicit yet. On the one hand, if two ele- 
Integrating it over a time T, one finds 

ments i and j belong asymptotically to two different orbits 

b and c, their distances (™ should remain finite. On the 

\xi{t + T) - Xj{t + T)\= exp(TAij) \xi{t) - Xj{t)\ , (4) other hand, since the phase space is finite, the distances 

Ay = In (a(l - e)) + — X! ^'^ l^^(^) + ^j(*)l ' (^) cannot diverge. Thus, the orbits of the two clusters have 
t 

to fulfill the condition 

If two elements belong to the same cluster &, their dis- 
tance has to shrink to zero. In this case, Xi{t) « Xj{t) « 

Xb{t). Kaneko ^T] defined the splitting exponent to mea- \c = In (a(l - e)) + ^lirn^ - ^ In \Xi,{t) + Xc{t)\ = . 
sure the rate of convergence to the orbit {^^(t)}, (7) 

For periodic orbits, this condition is just a consequence 
of periodicity. Nevertheless, it allows to rationalize some 



A, = In (2a(l - e)) + lim i ^ In \X,{t)\ , (6) 



T^oo 1 —^ features of GCLM. We call A^c repulsion exponent, since 

and defined an orbit as transversely stable if it has A^ < its positive value would mean that the two orbits repel 

(see also Q). While in the definition of the splitting ex- each other, and define two orbits as orthogonal if their 

ponent the partition iVi • • • Njc does not enter, it is the repulsion exponent vanishes. A set of orbits is stable if 

distribution of the elements into clusters which decides all of the orbits are transversely stable and all pairs of 

whether a set of orbits is a global attractor or not. All of orbits are orthogonal. This condition does not depend on 

the stable periodic orbits of the single map have negative the partition of the N elements into the K, orbits, but 

splitting exponent for every positive e. The splitting expo- a precise partition is needed so that the set of orbits is 

nent can be positive for non-entrained elements (forming invariant under the global dynamics. 
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Note that, for e = 0, the set of aU periodic orbits, stable 
and unstable, fulfill the orthogonality condition, which is 
just equivalent to periodicity. Thus, if the dynamics of the 
single map is periodic with period P, the P orbits obtained 
from the different phases of the stable periodic orbit con- 
stitute a stable set. As e increases, the transversal stability 
condition becomes easier to fulfill even for lower periodic- 
ity orbits, but the orthogonality condition becomes more 
difficult (notice that the larger the number of clusters, the 
more demanding this condition) . Thus at some point only 
partitions with less than P clusters can be found. This is 
probably a reason why only very small numbers of clus- 
ters are typically observed in simulations. Additionally, 
since the average value J2t ^^ \-^b{t) + Xc{t)\/T cannot be 
larger than log(2), no stable two cluster system can exist 
for e > 1 — l/2a. This is only an upper bound, since the ac- 
tual value of e where multiple cluster attractors disappear 
is much smaller. 

In Fig. we show the two splitting exponents as well 
as the repulsion exponent for a two-cluster system with 
parameters a — 1.3 and e = 0.15. There is a continuous 
spectrum of partitions allowed, and for all of them the 
clusters move along period- two orbits, which are periodic 
orbits of the two variable dynamical system 

Xi(f + 1) = 1 - a ((1 - e(l - p))X^,it) + eil - p)Xi{t)) 
X2it+l) = l-a{{l- ep)Xl{t) + epXlit)) , (8) 

where p = Ni /N is the fraction of elements in the largest 
cluster and is kept fixed during the dynamics. Notice that 
the repulsion exponent A12 equals zero up to very high pre- 
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Fig. 2. Possible two-cluster partitions for a system with a — 
1.3 and e — 0.15. For different values of the relative size p of 
the largest cluster the two clusters move along two period-two 
orbits represented as empty and solid circles, respectively. The 
fraction of initial conditions converging to two-cluster orbits is 
shown as a dashed line. The remaining trajectories synchronize 
completely to the stable period-four orbit. The two transversal 
exponents Ai and A2 as well as the repulsion exponent A12 are 
shown. Notice that Ai attains a minimum at p ~ 0.52. Empty 
circles represent the first period two orbit, full circles represent 
the second one. 

cision for all values of p. For these parameter values also 
the completely synchronized state is stable, and moves 
along a period-four attractor whose attraction basin cov- 
ers roughly 70% of phase space for the values of p where 
the two clusters are stable. At large p, the transversal 
exponent A2 approaches zero, and the two-orbit system 
becomes unstable while the synchronized attractor covers 
100% of phase space. 
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2.2 Replicas and overlaps totic attractors coincide. Even in this case, however, the 

above formula can take different values depending on the 

In [Q, two of us have introduced the concept of replica 

relative phase of the two orbits and on the permutation 

in GCLM and a measure of similarity between them, the 

of elements. To avoid this, the second orbit is shifted by 

overlap q"^ . A replica a is an orbit {xf{t)} of the whole 

a number time steps with respect to the first one in order 

system, and different replicas are obtained from different 

to maximize the overlap. This procedure is similar to that 

realizations of the initial conditions. Thus, the overlap qaf} 

used in spin glasses with rotational symmetry [|l5| . Finally, 

is a random variable dependent on the sets of initial condi- 

the degeneracy due to the arbitrary initial labelling can 

tions {xf{0)} and {xf (0)}. We investigate its distribution 

be avoided through a proper reordering of the elements 

for specific ensembles of initial conditions, keeping the pa- 

once the stable attractor has been reached: Elements in 

rameters e, a, and N fixed. 

the largest cluster are assigned labels from 1 to A^i, those 

In order to compute the overlap, we transform the or- 

in the second largest from A^i + 1 to A^i + A^2, and so on. 

bits into binary sequences, assigning a binary number ai{t) 

Using the previous definition, the overlap q returns a 

to each element i at each time step t such that cr'j'it) = 1 

finite positive value for clusters of periodic orbits and of 

ifxf{t) > X* , and erf (i) = — 1 otherwise, with 

chaotic ones in the two-band chaos, due to the regular 

alternation of plus one and minus one values in this re- 

X* = '^ + f+^ (9) 

2a gion. Orbits with one-band chaos change the sign of the 

the fixed point of a single logistic map. sequence a{t) in an uncorrelated fashion, implying that 

Finally, the overlap q°'^ is defined as in the limit T ^ oo their overlap with other orbits tends 

to zero. Our numerical simulations indeed show that this 

, to+T N 

q'^0 — y^ y^cr"(t)cr^(t) . (10) i^ the case. Thus, the definition of the overlap becomes 

problematic for orbits with one-band chaos. 



This quantity is computed after a large enough transient 



time to has elapsed, so that the two trajectories reach 
their asymptotic attractors, and averaged over the mini- 
mal common multiple of the two periods or, in case the 
asymptotic dynamics is not periodic, over a very large 
simulation time T ~ 10^~^, depending on the underlying 
dynamics. The overlap takes values between -1 and 1. In 
order for it to be a meaningful measure of similarity, the 
value 1 should be returned if and only if the two asymp- 



2.3 Attraction basin weights 

The overlap distribution gives information on the distribu- 
tion of attraction basin weights for the particular ensemble 
of initial conditions chosen. In fact, we can write it as 



F((Z)-^PF<,I^^%-g„0), (11) 

Q/3 
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where a and /3 label all possible global attractors, Qap is coherent attractor. It turns out that multiple-cluster at- 

their overlap, and Wa and W^ are their attraction basin tractors have vanishing attraction basins in the therniody- 

weights, i.e. the fraction of initial conditions which con- namic limit, so that the only nonvanishing contribution to 

verge to the attractors a and /3 respectively, and whose Y comes from the completely synchronized attractor, and 

sum is normalized to unity. The overlap distribution con- we can approximate Y ~ T^cs- Thus we can also study 

tains a delta distribution at q — 1, obtained for a = P, through the parameter y the transition between complete 

whose size is equal to the average attraction basin weight: and partial synchronization. 

Y^Y.^a- (12) 

a 3 Fixed-field ensemble 

This parameter expresses the probability that two ran- 
domly chosen initial conditions converge to the same at- To compute overlap distributions and other statistical prop- 
tractor and is able to distinguish between different situ- erties of GCLM, a set of replicas corresponding to differ- 
ations. If Y is equal to or tends to one in the thermody- ent initial conditions should be taken. Ideally, all initial 
namic limit, it means that there is only one relevant at- conditions should be present in the set. In an actual com- 
tractor which attracts in this hmit all of the phase space putation, this is never possible. Instead, a large ensemble 
of the system. If Y tends to zero in the thermodynamic of initial conditions is randomly generated. One expects 
limit, it means that the system has in this limit a diverg- that averaging over this ensemble would be equivalent to 
ing number of different attractors and none of them is the ideal averaging over "all" initial conditions. However, 
dominant. The situation in between, when Y is finite but this would only be true if the employed random set is rel- 
smaller than one, means that there is a finite number of atively uniformly sampling the full space of initial condi- 
relevant attractors. tions. Some random sets of initial conditions may be miss- 
For a fixed value of a and increasing e, the collective be- ing this property. In previous numerical investigations [ |lO[ 
haviour of GCLM changes from turbulent to glassy (mul- O] of dynamical glasses, to generate a set of initial con- 
tiple clusters) to finally fall into complete synchronization, ditions the coordinate 2:^(0) of each map J = 1, ..., N was 
This last transition can be characterized through different chosen independently and with a uniform probability den- 
parameters. Kaneko proposed to use the average cluster sity from the interval (—1, 1). It was tacitly assumed that 
number |l| , which grows continuously from a finite value this procedure would yield uniform sampling of initial con- 
in the dynamical glass phase to unity in the synchronous ditions. However, the initial conditions generated in this 
phase (for fixed N). An alternative measure can be the way become increasingly similar in the thermodynamic 
fraction Wqs of initial conditions converging to the single limit N —i- oo. 



8 S.C. Manmbia, U. Bastolla, and A.S. Mikhailov: Glassy behaviour of GCLM 

As follows from (m and ^ , the dynamics of GCLM is Below in this section we discuss glass properties of 

described by the equations GCLM for evolutions starting from a fixed field ensem- 

ble (to our knowledge, this is the only way in which the 
x^{t + 1) = 1 - a [(1 - e)x^^{t) + em{t)] (13) 

initial conditions have been so far modelled in the litera- 



■^here ture) . As we shall see, this ensemble leads to replica sym- 

Tn{t) = — 2_^ ^i \t) (14) metry, since all approached attractors are then identical 



N 

i=l 

is the synchronizing field that acts on a given element z ^P *« ""1^11 variations which vanish in the thermodynamic 
and is collectively produced by the whole system. miit. 

Let us consider statistical properties of the initial syn- 

, n 1 , /, n\ ■ ^-u ^■ -I- f ^ AT 3.1 Transition to complete synchronization 

chromzmg neid m — m[t — U) m the umit of large iV, r j 

when the coordinates Xi{0) of each map i — 1, ...,7V are 
chosen independently and with a uniform probability den- 
sity from the interval (—1, 1). Because this field represents 
then a sum of a large number of independent random vari- 
ables, it should obey for N ~* oo a, Gaussian probability 
distribution 



p(m) = , exp 

V2w 



{m — m) 



(15) 



2cr 

where m is the mean value of the field m and a is its 
mean-square statistical variation. Using (|l4|), we obtain 

N -, N 



At sufficiently high coupling strength e, GCLM become 
completely synchronized. We characterize the transition 
to complete synchronization through the probability Y = 
P{q = 1) (see Eq. (p^)) that two randomly chosen initial 
conditions from a fixed-field ensemble fall into the same 
attractor, see Fig. 3. Our analysis is performed with a 
value of the logistic parameter a = 1.3 where the single 
logistic map is periodic with period four. The coexistence 
of different attractors and their dynamics for this param- 
eter choice have been previously considered in M. Near 



1 1 fi 1 

77j — y^(a;^fO)') == Y^ / x^dx = - (16) the transition to complete synchronization the final stable 

^/1^\ 2\ J/ 2N'^J_i 3 ^ ^ 

i=l 1=1"'^ 

attractors consist of two clusters following the dynamics 
and 

of period two. 

1 ^ A 

(7 = ((m — m)^) = -j^ /^((^t(0)) ~ (2^i(0))^) — .- ., It is interesting how this transition takes place in the 

i—l 

(17) thermodynamic limit. Each curve in Fig. 4 represents Y 

Thus, in the thermodynamic limit N ^ oo the initial for a given value of e as a function of N . Even if the 

synchronizing field approaches a constant value, indepen- completely synchronized state is stable for e > ei(a), the 

dent of the realization. For large A^'s it shows fiuctuations smallest value at which the synchronized state is transver- 

of order I/a/A. Such a set of initial conditions shall be sally stable, this state is never observed until e reaches a 

called a fixed- field ensemble below. larger value e > £2(0) |Q. The transition is discontinuous 
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Fig. 3. Transition to complete synchronization. When cou- 
pling e is increased, the attraction basin of the completely syn- 
chronous state grows until all initial conditions end there. 
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Fig. 4. Average weight of an attractor basin Y = P(l) as a 
function of the system size for a = 1.3. 

(first order) and takes place at e ~ 0.165 for a — 1.3. For 
smaller e the system is in a phase where many different 
attractors coexist. All of them are two-cluster attractors. 
Since the average attraction basin weight Y vanishes in 
the thermodynamic limit, the number of different attrac- 
tors diverges for A^ — > oo. The unbounded increase of the 
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Number of different initial conditions 

Fig. 5. The number of different attractors reached by the sys- 
tem increases as a power law of the number of random initial 
conditions used. Increasing the latter number is equivalent to 
exploring in higher detail the phase space. From top to bottom, 
iV = 16384, 4096, 1024, and 256. 

number of different attractors for iV — > oo was indeed 
one of the first indications that GCLM might represent a 
glass- like system Q. 

We have further examined how the number of attrac- 
tors Ai visited by the system grows as the number I of 
different initial conditions used increases. Our results for 
e = 0.15 and a — 1.3 are displayed in Fig. 5. We observe 
an approximate power-law dependence Al ex X'' , with an 
exponent rj dependent on the system size A^. For 2 — > oo, 
Ai should saturate at a finite value. Although the number 
of different attractors grows fast, the bending at large T 
for the largest size reflects the existence of an asymptotic 
value A^oo(A^)- 
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Fig. 6. Data collapse of the normalized overlap distribution 
under increasing system size; system parameters e = 0.15 and 
a = 1.3. 

3.2 Distributions of overlaps 

To quantify the similarity between different attractors, we 
have calculated the overlap distributions for the same pa- 
rameters, e ~ 0.15 and a — 1.3, and different system sizes 
N . As seen in Fig. 6, such distribution P{q) approaches a 
delta- function in the thermodynamic limit A^ — > oo. The 
width of the main peak in P{q) goes to zero as 1/vN. 
Some finite-size effects can be observed in the bump at the 
smallest size represented, and in the peaks which appear 
intermittently for relatively large values of N, showing the 
"locking" of the system close to prefered partitions, before 
reaching the asymptotic behaviour. 

The finite-size behaviour can also be more compli- 
cated. Fig. 7 shows the overlap distributions obtained at 
the same control parameter a ^ 1.3 for four different val- 
ues of e and systems of size N = 256. For small e (Fig. 7a), 
there is a large number of partititions close to the attrac- 



tor with the largest basin, Ni — 149, ^^2 = 107. A second 
group of attractors corresponds to three-cluster families 
close to A^i == 115, N2 == 85, N3 = 55. The attractors 
within each group are similar and their mutual overlaps 
are close to unity, explaining the large weight of P{q) at 
q ~ 1. The overlaps between these two groups give a sec- 
ond contribution around q = 0.7. The continuous line in 
Fig. 7a shows the total distribution, the dashed line corre- 
sponds to the attractors with the same number of clusters, 
and the dotted line represents the contribution from the 
overlaps between two- and three-cluster attractors. The 
part of the distribution close to q = I results from three- 
cluster attractors where the third cluster has only a few 
elements. As the coupling strength grows, the three-cluster 
attractors become less and less frequent, and two-cluster 
attractors dominate (Fig. 7b, c). For large enough coupling 
(an example is e == 0.17 in Fig. 7d), the completely syn- 
chronous state appears and starts to occupy an increas- 
ingly large fraction of the phase space. Its self-overlap is 
unity, while its overlap with the remaining two-cluster at- 
tractors is small (three-cluster attractors are no longer 
present). The overlap between one- and two-cluster at- 
tractors explains the large contribution at small values 
of q observed in Fig. 7d. If e increases further, the com- 
pletely synchronous state attracts more and more initial 
conditions, P{q) tends to a delta-function at g = 1, and 
the contribution at small q disappears. 

In the chaotic domain of a single logistic map, for 
a > Coo, a similar behavior but with strong finite-size 
effects is observed. In a previous publication |12|, we have 
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l/\/N . The origin of such fluctuations lies in the 1/vN 
variations of the synchronizing field. 
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Fig. 7. Overlap distributions for A'^ — 256 and a — 1.3. (a) 
e = 0.05, (b) e = 0.1, (c) e = 0.14, and (d) e = 0.17. 



studied the parameters e = 0.1 and a = 1.55. This point 
had been also analysed in fin] for its glassy properties. The 
overlap distribution is broad here and its shape keeps al- 
most unchanged with the system size until N ~ 4000. But 
when N increases further, the most of the distribution's 
weight is shifted towards q — 1, indicating that the attrac- 
tors reached by the system indeed become very similar. 
For e = 0.3 and a — 1.9, we have observed that for sizes 
up to A^ ~ 3000, the overlap distribution P{q) remains 
almost constant (see Fig. 4 in ll2J). If we apply rescaling 
similar to Fig. 6, only the two largest sizes {N = 2048 
and N ~ 8192) seem to follow the expected asymptotic 
behaviour and collapse. 

Our analysis shows that, for N large enough, GCLM 
tend to a prefered cluster partition. It can be said that the 
same macroscopic state is always found in the thermody- 
namic limit, apart from "thermal fluctuations" of order 



3.3 Distributions of cluster sizes 

Information similar to the overlap distribution is contained 
in the distribution Q{pk) of cluster sizes pfc — Nk/N, i.e. of 
the fractions of elements belonging to a cluster (obviously, 
X]fc-Pfc ~ !)■ For large sets of initial conditions (varying 
between lO'' and 2 x 10^ realizations) in the fixed-field en- 
semble, we have computed the values of pk for all stable 
partitions and thus obtained the distributions Q{pk) of 
cluster sizes. 

An example of such distributions for e = 0.15 and 
a = 1.3 and different system sizes is shown in Fig. 8a. We 
see that in this case the asymptotic attractors are always 
formed by two clusters of unequal size. Their dynamics 
is periodic with period two. As iV — ^ oo, the distribu- 
tion Q{pk) shrinks in width around the two prefered sizes, 
pi ~ 0.372 and p2 = \ — pi- For large enough N the two 
peaks approach a Gaussian distribution, and its width de- 
creases proportional to l/vW , as shown in Fig. 8b. 

A similar behaviour was observed for other parameter 
values. Generally, for sufficiently large TV the distribution 
of the sizes of the fcth cluster is a Gaussian centered at a 
prefered value p^. We show two more examples. In Fig. 9a 
(e — 0.3 and a — 1.9), the system tends to attractors with 
two clusters of almost equal sizes, though occasionally also 
three-cluster attractors are observed. Fig. 9b (e = 0.1 and 
a — 1.55) gives an example where attractors with four 
clusters of different sizes are prefered. We show the total 
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Fig. 8. (a) Distribution of cluster sizes for increasing system 
size; parameters e = 0.15 and a — 1.3. (b) Data collapse of the 
size distribution for the largest cluster, with A'^i ~ 0.628A'^. 

size distribution Q{pk) together with the size distributions 
for clusters of rank one to four. Here, the prefered partition 
is close to pi = 0.307, p2 = 0.242, ps = 0.2295, and p4 = 
0.2215. 

Thus, we have investigated numerically the statistical 
properties of GCLM in the glass phase starting from initial 
conditions in the fixed-field ensemble. The investigations 
show that in the thermodynamic limit A^ — > oo this system 
has a great number of different attractors, increasing as a 



power law of system size N. However, all these attractors 
are very similar. Namely, the differences in their statistical 
properties, such as the cluster sizes, are proportional to 
I/VN and thus vanish in the limit N -^ oo. This explains 
why replica-symmetry is recovered in the thermodynamic 
limit, when only the evolutions initiating from a fixed-field 
ensemble are considered. 

4 Random-field ensemble 

4.1 The role of initial conditions 

For given parameter values a and e, many different or- 
bits corresponding to a continuous spectrum of two-cluster 
partitions and to complete synchronization are stable (see 
Fig. g). Yet, only one partition is chosen in the fixed-field 
ensemble, up to variations of order 1/yN. The selection 
of this prefered partition cannot be explained by a higher 
stability of its orbit. For instance. Fig. || indicates that for 
a = 1.3 and e = 0.15 the transversal stability is strongly 
increased at the cluster size pi ~ 0.522. However, the se- 
lected partition in the fixed-field ensemble has in this case 
the cluster size pi ~ 0.628 whose stability is much weaker. 
For parameters in the chaotic domain of a single logis- 
tic map, the situation is similar. Stable attractors with 
chaotic dynamics exist here, but the system often shows 
a preference for the periodic ones. 

To examine in more detail the role of initial conditions, 
we use a slightly generalized version of the fixed-field en- 
semble. Namely, we assume now that the initial coordi- 
nates Xi{0) of all maps are independently and uniformly 
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limit A^ ^ oo it is again given by a Gaussian distribution 
with mean value m == ^/3 and mean-square dispersion 
(7 = (4/45) (^/A^). Hence, at time i = 1 the coordinates 
Xi{l) of the maps are given by 



x,(l) = 1 - a [(1 - e)a;,2(0) + em 



l-a 




(l-e)x?(0) + 6| + O(7V-i/2) 



(18) 
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Fig. 9. (a) Distribution of cluster sizes for different system 
sizes, e = 0.3, and a = 1.9. For A*' large enough, the sys- 
tem chooses a partition formed by two clusters of similar sizes, 
(b) Distribution Q{pk) for N = 8192 and parameters e — 0.1 
and a — 1.55. The whole distribution Q{pk) and individual 
distributions for the largest, second, third, and fourth largest 
clusters are displayed. 

distributed between — ^ and ^ (note that the fixed-field 
ensemble corresponds then to the choice ^ = 1). Cal- 
culating again the statistical distribution of initial syn- 
chronizing fields TO, we find that in the thermodynamic 



Since 2:^(0) is uniformly distributed, xf{0) is distributed 
with density l/2a; which diverges at a;^(0) = 0, so that 
the distribution of Xi{l) has a pronounced maximum at 
X = I — ae^/3 . This initial bias drives the system to- 
wards the attractor closest to the most probable value of 
Xi{l). This is shown in Fig. Wu, where we represent, as 
a function of the fixed ensemble parameter ^, the value 
1 — ae^/3 where the maximum of 2:^(1) is expected, the 
actually observed maxima of a;i(l), and the coordinate on 
the prefered asymptotic orbit. Varying ^, the bias in the 
initial value 2:^(1) changes and drives the system to dif- 
ferent asymptotic orbits (full circles in Fig. HG), in turn 
corresponding to different partitions of the elements. We 
thus find p = 0.632 at ^ = 1 and p = 0.717 at C = 0.944. 
Partitions with larger p cannot lead to two transversely 
stable orbits, thus for ^ < 0.944 the completely synchro- 
nized attractor is always reached. 

It is interesting, in this framework, to look at the dy- 
namics of synchronization for a = 1.3, e = 0.15 and vary- 
ing ^ (see Fig. Ill|). Different examples always show the 
same pattern: after the first time step, the most popu- 
lated region of phase space coincides with the maximum 
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Fig. 10. For different values of the initial field, the figure shows 
the most likely value of Xi{l) (empty circles) and the coordinate 
of the prefered attractor closest to this value (full circles) . For 
£, < 0.944 the prefered attractor is completely synchronized. 

in the distribution of Xi{l). After very few time steps, 
the two most populated "clusters" start to oscillate close 
to a stable attractor made of two periodic orbits of pe- 
riod two, but most elements are not synchronized yet and 
the partition is quite different from the final one. At the 
same time as oscillations around the periodic orbits are 
dumped, more and more elements join the two clusters, 
until the partition which stabilizes the periodic orbits is 
reached. Thus the system first chooses the orbits and only 
afterwards partitions which would stabilize them. 

A similar route is observed even when the system tends 
to the completely synchronous period- four orbit. First the 
system approaches two period-two orbits which are very 
close one to each other and starts to partition on them. At 
some point, the smaller cluster is attracted by the larger 
one and disappears. For some parameter values the period- 
two orbit remains metastable for quite a long time, until it 



Fig. 11. Dynamics of two-cluster synchronization for a = 1.3, 
e = 0.15 and A*' — 1000. The solid circles represent the coordi- 
nates of the largest cluster, empty circles those of the second 
larger cluster, the solid and dashed lines represent their sizes 
respectively. 

splits into a period- four orbit through a kind of dynamical 
bifurcation (see Fig. |l2|). Even for values of a in the chaotic 
phase we observed synchronization first through attrac- 
tion towards prefered period-two orbits and then through 
a bifurcation (see Fig. [l3| ). 

Summaryzing the findings of this section, we can say 
that the initial synchronization field strongly biases the 
elements towards a prefered region of phase space, lead- 
ing them to periodic orbits which can be either stable 
(and indeed stabilized through the appropriate partition 
of the system) or metastable (and eventually transformed 
to completely synchronous attractors). 
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bifurcation, it reaches the stable period-four orbit. 

4.2 Replica-symmetry breakdown in the random-field 
ensemble 






As shown in the previous section, by varying the parame- 
ter f we can drive the system to macroscopically different 
attractors. Thus an ensemble of initial conditions, where 
^ is randomly chosen for each initial state, is expected to 
lead to very different attractors. We define a random-field 
ensemble as a random set of initial conditions which is 

generated in the following way: For each realization, we Fig. 14. Distributions of overlaps for a = 1.3 and e = 0.15 
first choose at random the parameter i form the interval ^^^ ^^e random-field ensemble for different system sizes iV. The 
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(0,1). Then the initial states Xi(0) of all individual maps 
in the system are independently drawn from the interval 



value of P(l) is not represented. It accumulates around 50% 
of the total weight of P(g). 
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Fig. 15. Distributions of overlaps for a = 1.55 and e = 0.1 in 
the random-field ensemble for different values of TV. For these 
parameters the single map is chaotic. The value of P(l) — 
for all N. 

When such random-field ensembles of initial condi- 
tions are used, overlap distributions do not shrink to a 
delta-function peak, but remain continuous in the ther- 
modynamic limit. This is shown in Fig. 11^, which dis- 
plays overlap distributions in the random-field ensemble 
for a — 1.3 and e — 0.15 and different system sizes. In this 
case, the weight of the completely synchronized attractor, 
Y ~ P{q — 1) does not vanish in the glassy phase. We ex- 
pect that the transition to complete synchronization is in 
this case second-order like: Y tends continuosly to unity 
as e approaches the critical coupling at which only syn- 
chronized orbits are stable. 

We present also the distributions of overlaps for pa- 
rameters in the chaotic domain of the single map, e = 0.1 
and a = 1.55, and three different values of the system 
size N. There is again a qualitative difference between the 



fixed- field and the random- field ensemble. While in the 
former case the function P{q) showed a systematic loss of 
structure for increasing N (see Fig. 1 in [f2|), in the latter 
situation it remains remarkably invariant with the growth 
of the system size. 

Thus, we see that for the random-field ensembles of ini- 
tial conditions the overlap distribution becomes indepen- 
dent of the system syze in thermodynamic limit N —^ oo . 
The asymptotic overlap distribution is formed by a delta- 
peak at g = 1 plus a broad, smooth part extending to low 
values of the overlap q. The presence of such continuous 
distribution is an indication of replica-symmetry break- 
ing. The breakdown of replica symmetry means that, for 
each orbit of GCLM, one can find orbits of gradually vary- 
ing degrees of similarity within a large ensemble of orbits 
generated by randomly chosen initial conditions. 

4.3 Ultrametricity 

Generally, the ultrametric distance d{A, B) between two 
elements A and B in a hierarchy is defined as the num- 
ber of steps one should go up in the hierarchy to find a 
common ancestor of two elements A and B. If any three 
elements A, B and C belong to a hierarchy, the inequality 
d{A,C) < ma.x{d{A,B),d{B,C)} should hold. As a con- 
sequence, the two maximal distances between elements in 
any triad must always be equal. If overlaps q"^ between 
any two replicas a and (3 are uniquely determined by the 
ultrametric distance d{a, (3) between the respective states, 
the overlaps between any three replicas a, (3 and 7 must 
satisfy the relationship 
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Fig. 16. Distribution H{Aq) of distances between the two 
smaller overlaps our of a triad. In the thermodynamic limit 
this function has a non-vanishing width. Same parameters as 
in Fig. 15. 



^''>min{g"^,g"''}, 



(19) 



implying that the two minimal overlaps in any triad of 
replicas are always equal |]l3J . 

To check the presence of ultrametricity, we have to 
consider triads of replicas a, (3, and 7 and calculate the 
three overlaps that can be defined by combining them. If 
the two minimal overlaps in any triad are always equal, 
the ultrametricity is present. Formally, this amounts to 



requiring that the relationship (19) always holds. That 
condition can be numerically tested by generating triads 
of replicas and computing the distribution H{Aq) over 
the differences between the two minimal overlaps, Aq = 
l^a/3 _ ^a7|^ jf H{Aq) -> 5{Aq) in the limit TV ^ 00, then 
the system is ultrametric. 



Previously, such calculations have been performed for 
the fixed-field ensemble [|2). In this case case the distri- 
bution H{Aq) approaches a delta-function 5{Aq) in the 
limit of large system sizes N . However, as becomes clear 
from the analysis of overlap distributions in the present 
study, this behaviour simply refiects the vanishing diver- 
sity of system attractors for the fixed-field ensemble in the 
thermodynamic limit. 

We have now repeated such calculations for the random- 
field ensemble. Distributions H{Aq) of distances between 
the two smaller overlaps in randomly generated triads of 
replicas for two different system sizes are shown in Fig. 16. 
We see that the distributions are broad and almost do not 
depend on the system size. Thus, GCLM do not display 
ultrametric properties. 

Though replica-symmetry breaking is a necessary con- 
dition for nontrivial overlap distributions, it does irot im- 
ply exact ultrametricity, which is a much more demand- 
ing condition. Possible deviations from exact ultrametric- 
ity have been discussed for spin glasses [|l3|. Parisi and 
Ricci-Tersenghi [Oj have shown that exact ultrametricity 
can only hold under the conditions of stochastic stability 
(i.e. that each replica is in a certain sense equivalent to 
the others) and of separability (i.e. that all the mutual 
information about a pair of equilibrium configurations is 
already encoded in their overlap). Our nunrerical analy- 
sis of GCLM shows that for the random-field ensemble in 
the thcrmodynamical limit this system is characterized by 
replica-symmetry breaking, but exact ultrametricity is ab- 
sent. Note that though exact ultrametricity, which would 
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have corresponded to the appearance of deha-function coupUngs larger than the critical one, the system synchro- 
peak at Aq — 0, is not observed, the distributions H{Aq) nizes completely for nearly all initial conditions in the 
in Fig. 16 have a broad maximum at Aq = 0. This indi- thermodynamic limit, and the average attraction basin 
cates that some weaker form of organization may still be weight tends to unity. This situation is reminiscent to the 
present here. analogous transition in attraction basin weights observed 

for random boolean networks |16[ and for asymmetric neu- 

5 Discussion and conclusions ral networks 0. In both cases, the average attraction 



We have examined asymptotic glass properties of glob- 
ally coupled logistic maps in the thermodynamic limit for 
two different random ensembles of initial conditions. In 
the fixed-field ensemble the initial value of the synchro- 
nizing field becomes identical up to variations of order 
1/ yN that vanish in the thermodynamic limit. Therefore, 
all attractors reached by the system become increasingly 
similar for N —> oo. Dynamically, the bias due to the ini- 
tial field drives the system towards the prefered attractor 
and then the elements partition in such a way to stabilize 
the prefered attractor. The overlap distribution tends to 
a delta- function peak at g = 1, i.e., even when a diverging 
number of attractors is present, they are all macroscop- 
ically identical. Hence, replica symmetry is recoverd for 
the fixed-field ensemble in the thermodynamic limit. 

We have also found that, in the fixed-field ensemble, 
the system undergoes a special phase transition when e 
overcomes a critical value. For smaller coupling, the sys- 
tem is partitioned into a small number (close to the tran- 
sition, usually two) of periodic orbits. Though all of the 
attractors reached for different initial conditions are very 
similar, their number diverges in the thermodynamic limit, 
and their average attraction basin weight goes to zero. For 



basin goes discountinuously (in the thermodynamic limit) 
from the value zero, when the system is in the "ordered 
phase" , to a finite value related to the average attraction 
basin of random maps pq | in the chaotic phase. Though in 
GCLM the finite attraction basin weight is a consequence 
of complete synchronization, the formal analogy between 
this system and dynamical systems with quenched disor- 
der is very suggestive. 

The asymptotic behaviour of GCLM in the thermody- 
namic limit is essentially different when the random-field 
ensemble of initial conditions is chosen. Because initial 
synchronizing fields retain in this case macroscopic fluc- 
tuations even for iV — > cx), a broad range of attractors may 
still be reached. In the random-field ensemble, the tran- 
sition to complete synchronization (with the average at- 
traction basin weight F = 1) is expected to be continuous, 
more in analogy with equilibrium mean-field spin glasses. 
Examination of the overlap distributions has revealed that 
replica-symmetry breaking persists in the thermodynamic 
limit. Thus, GCLM reach the status of a dynamical coun- 
tepart to mean-field spin glasses. 
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